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The increasing demand of energy and development of distributed generation 
systems have led to great progress in medium voltage high-power inverters. 
On the other side, due to more intransigent regulations in the area of power 
quality and appearance of new grid codes, which limit the harmonics’ 
amounts in ac grids, special attention has to be made to restrict the 
harmonics. The recently proposed selective harmonic mitigation pulse width 
modulation method reduces the values of specific harmonics to meet the grid 
codes, instead of making them zero. In this way, there are much more 
solutions for mitigating the same harmonics which in the Selective Harmonic 
Eliminated Pulse-Width Modulation (SHEPWM), are completely eliminated. 
The SHMPWM was first introduced in, where a waveform having 750 Hz 
switching frequency was employed to mitigate harmonics up to 49th to 
satisfy the grid codes EN 50160 and CIGRE WG 36-05 without using any 
filtering system. In this contribution a solution to the parallel resonance 
problem that can be present in practical applications of shunt Active Power 
Eilter (APE) compensation is proposed. The proposed solution involves 
turning the shunt APE scheme of compensation into a Shunt Hybrid Eilter 
(SHF) configuration. A Linear Quadratic Regulator (LQR)-based switching 
controller was specifically designed for this hybrid scheme of compensation, 
maintaining stringent performance requirements on the tracking of filtering 
currents and the draining of the harmonic ripple currents. 
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1. INTRODUCTION 

With the continuous increase in the incorporation of nonlinear loads in the power network, the 
harmonic distortion of voltage and current waveforms, among other power quality problems, has become a 
great concern. This paper focuses on the use of shunt filtering techniques such as shunt Active Power Eilter 
(APE) and shunt hybrid filter (SHE), to mitigate harmonic distortion and to compensate reactive power. In 
conventional inverters which use sinusoidal pulse width modulation, there are heavy losses due to presence 
of harmonics [1]. Upcoming smart grid paradigm enabling maximum power throughputs and near- 
instantaneous control of voltages and currents in all links of the power system chain [2]. The series capacitor 
provides the majority of the power pulsation decoupling through a wide voltage swing, and the buffer 
converter only needs to process a small fraction of the total power of the entire architecture, allowing a very 
small active circuit volume and very high system efficiency [3]. Inverter effectiveness is evaluated by 
comparing the performance based on SHE PWM and based on the conventional constant switching frequency 
PWM technique [4]. By installing DG in the transmission system, voltage stability and voltage profile can be 
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improved, while losses can be minimized [5]. In particular, the shunt APF has been demonstrated to be an 
appropriated tool for the mitigation of harmonic currents and reactive power compensation [6]-[15]. The 
shunt APF can be modeled as a controlled current source that supplies a compensation current in parallel 
with the nonlinear load. The main components of a shunt APF are: a reference current generator, a switching 
current controller based on a pulse-width modulation (PWM) technique, a power electronic converter 
(usually a voltage-source-type PWM inverter), and a dc voltage controller (PI controller). Usually, the 
modulation technique used when modeling the shunt APF is the hysteresis-band current controller, due to its 
fast dynamic response and its easy implementation [16]-[18]. 

The path should be provided to drain this high frequency ripple current; this goal could be achieved 
with a LC ripple-filter [19] or a compensator-passive filter structure designed for this particular purpose [20]. 
The necessity of adequate development an application of techniques to solve potential problems of resonance 
and harmonic distortion in diverse electrical networks has been stressed in more recent contributions [21]- 
[23]. In [24] frequency scan and modal analysis methods are applied for harmonic distortion and system 
resonances identification in wave power plant applications, so that appropriate correction measures can be 
timely undertaken. A distributed series/hybrid-shunt compensation technique is proposed in [25] for 
harmonic mitigation in industriaFcommercial facilities. The main aim of this paper is to provide a strong and 
high quality scheme of compensation, e.g., shunt filtering, when applied to non-stiff systems. 

The proposed control is good enough to provide the following goals: (1) To provide a sound 
solution to the parallel resonance problem that can be present in shunt APF compensation, (2) To allow the 
correct shunt filtering compensation, (3) To avoid the propagation of the harmonic ripple current into the 
electric network, (4) To reduce, as much as possible, the KVA rating of the APF. 


2. RESEARCH METHOD 

2.1. Shunt hybrid filter - parallel resonance phenomenon 

A major drawback of active power filters is their high rating, e.g. up to 80% of the nonlinear load in 
some practical applications [26]. The problems emerge when the APF becomes a costly solution for 
harmonic current mitigation and reactive power compensation. Due to economic concerns and problems 
associated with insulation, as well as series or parallel connection of switches [27], different structures of 
hybrid filter have evolved as a cost effective solution for the compensation of nonlinear loads [28]-[29]. The 
hybrid filter topology presented in this paper consists of a shunt APF, and a shunt capacitor bank acting as a 
reactive power compensator. For this compensation scheme, the shunt APF placed to perform harmonic 
elimination and reactive power compensation is not included in its compensation mechanism. 


2.2. Shunt hybrid filter control 

The LQR method is a powerful technique for designing controllers for complex systems that have 
stringent performance requirements. The LQR method seeks to find the optimal controller that minimizes a 
given cost function. This cost function is parameterized by two matrices, Q and R, that weight the state 
vector and the system input, respectively. A LQR switching controller is designed to be able to generate the 
control signals of the voltage source inverter of the shunt APF under the SHF compensation scheme. It can 
be observed how the reference current generator and the LQR-based switching current controller are 
incorporated together to complete the proposed control strategy for the SHF compensation scheme. The flow 
chart of the LQR calculations to obtain the feedback gain matrix K. Given a set of weighting matrices, Q and 
R, the process in the calculation of matrix K becomes simple. Most of the times, the weighting matrix Q 
reflects the importance of the states to be controlled and the selection of an appropriate value of R has to be 
chosen to provide the control effort of the system and to obtain a finite time convergence of the regulator. 
Once the LQR algorithm is solved, and the feedback matrices are obtained.The reference current generator to 
be used by the APF is implemented according with [5], and is given by. 


ifa(t) 

4(t) 

ifc(t) 


ila(t) - 


iib(t) - 


iic(t) - 



(1) 

2(PT+Ploss)sin((^^ + C(,^^-^) 

(2) 

2(PT+Ploss)sin((^t+C(,^^+^) 

(3) 


where i^(t) is the reference filtering current of phase k, iik(t) is the line current of phase k, e.g. in this 
example it corresponds to the load current drawn by the combination of the linear unbalanced load current 
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and the rectifier load current; the subscript k represents the phases a, b or c, respectively. Pt is the total active 
power delivered to the load, Ut is the sum of the peak value of the source voltages uk, /al is the phase angle 
of the fundamental component of the load voltage at phase a, and x is the angular speed at fundamental 
frequency f of 60 Hz. The capacitor reactance and capacitance are calculated as 


Vll 

(4) 

Q3c|> 

1 

(5) 


where Xcf is the capacitive reactance, Vll is the line-to-line rms voltage, and Cf is the capacitance per phase. 


Parallel resonance phenomenon of Load voltages and source currents: to identify the possible 
resonance problem, an analysis using the equivalent circuit of the electric system for phase a is performed. 
The linear load impedance is also included in the analysis. The equivalent impedance of the circuit at any 
harm harmonic order h is given by. 


(l/(Rs+jha)Ls))+ (l/(Ra+jha)La))+ (l/(Hjp^)) 

The test system of Figure 1 was implemented in Matlab/Simulink and the simulation results are 
shown in Figure 3. Figure 3(a) illustrates the load voltages; once the shunt APF is switched on after one 
period. Discordant waveforms appear, and this behavior may be caused by a resonance condition produced 
by the harmonic interaction of the APF, the capacitor bank and the load. Consequently, the same situation is 
presented in the source current waveforms of Figure 3(b). 


3. LQR-SWITCHING CONTROLLER 

For the designed LQR-based switching controller the optimum control of the system is based on 
minimizing the performance index of the form [12], 

I = - Wref)'^Q(w - Wref) + U^Rujdt (7) 

where w is the state vector and w^f is a desired reference state vector, Q is a Hermitian matrix or a symmetric 
positive semi-defmite matrix, R is a symmetric positive definite matrix. For this particular case, R is scalar 
and sets a penalty on the maximum control action, m is the optimal control law and T is the transpose 
operator. The minimization of the cost function results in a feedback control of the form 

m = -K(w - Wref) (8) 

where the gain matrix K is obtained by the solution of the steadystate Riccati equation [7]. Once a suitable 
gain matrix K is selected, the switching control is obtained from [12],[15]: 

m = -hys(K(w - w^ef)) (9) 

where the hys function is defined by 

if K(w — Wref) > ^ then hys(K(w — Wref)) = 1 (10) 

else if K(w — Wref) < ^ then hys(K(w — Wref)) = —1 (11) 

This technique is called switching band tracking control in [12], and the selection of Di determines 
the switching frequency while tracking the reference. The equivalent circuit of the hybrid filter compensation 
scheme, where the combination of both linear and nonlinear loads is represented by a series RL equivalent 
load. The state vector x and the state-space equation are obtained from the analysis. 

( 12 ) 
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For convenience, the estate vector x can be transformed into a new state vector w, which is 
expressed in terms of the local variables of compensation and/or the control variables, e.g., the filtering 
current if, the current trough the capacitor ic, the terminal voltage of the capacitor uc and the load current il. 
From the relationship between the state vector x in (11) and the new state vector w is observed. The state 
transformation is finally made by using the transformation matrix V as shown in (13), 
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(15) 


For the control law represented by (8), it is assumed that we have full control over m and the 
switching control is based on a linear combination of multiple states [15]. For most realistic applications, the 
LQR problem is solved via a computational package, such as Matlab. The Matlab command IqrdK; Q; RF, 
with 14 1 Jt 2, solves the LQR problem and the associated Riccati equation. This command calculates the 
matrix of optimum feedback control K, so that the control feedback law (8) minimizes the performance index 
(7) subject to (14). Most of the times, the weighting matrix Q reflects the importance of the states to be 
controlled and the selection of an appropriate value of R has to be chosen to provide the control effort of the 
system and to obtain a finite time convergence of the regulator [12]. An alternative for the tuning decision of 
the weighting matrices is to analyze the step response of the system in order to jump-start the learning curve 
necessary to select Q and R and be able to apply the LQR design method to a practical problem. Once the 
LQR algorithm is solved, and the feedback matrices are obtained, the implementation would proceed with no 
further calculations of the LQR algorithm, unless the circuit and/or load conditions are changed. 


4. RESULTS AND ANALYSIS 

The Shunt Active Power Filter is implemented in MATLAB is shown in Figure 1 and the 
corresponding non linear Load is shown in Figure 2 and results are described. For instance, the weighting 
matrix to minimize the performance index of the LQR problem is selected as Q 14 diag QVi 25 10 3 0 F, 
where diag stands for diagonal matrix, and R 14 0:01. Although the LQR switching band current controller is 
claimed to be robust, and sometimes, only the calculation of the feedback matrix of a single-phase is 
performed and the results used for the other two phases, for this case study, we have decided to calculate the 
reduced feedback gain matrices for each phase a, b, c, according the flow chart. The resulting feedback 
matrices for each phase are given as, 

K, = [ 49.9468 9.3787 34.5637] 

Kb = [ 49.9460 9.3795 34.7046] 

K,= [49.9450 9.3805 34.7665] 

The results obtained from the simulation with Matlab/Simulink are shown in Figure 3. It can be 
observed that now, both the load voltages shown in Figure 3(a), and the source currents shown in Figure 3(b), 
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become balanced and free of harmonics once the shunt APE goes into operation at one cycle of simulation. 
The THD index for load voltages and source currents is null. In addition, the source currents are placed in 
phase with the load voltages, thus compensating the reactive power and achieving a unity load power factor. 
In similar way, from Eigure 3(a), Eigure 3(b) and Eigure 3(c), that there is not harmonic ripple contamination 
into the electric network. The high-frequency ripple current has been drained-out directly through the 
capacitor branch avoiding its propagation into the electric network. Most importantly, the parallel-resonance 
problem. This goal has been achieved despite that the frequency response of the electric system is the same. 
The shunt APE current injected at phase c is illustrated in Eigure 3(d). As a comparison exercise please 
observe the current injection in phase c of Eigure 3(d) and compare it to the current injection shown in Eigure 
3(d), which is the current injected when the compensation of reactive power was planned to be achieved by 
the shunt APE. It can be also noticed that the ripple current is now smaller in amplitude. Eor this case, the 
three-phase rams filtering current was 8.54 A, giving a KVA rating of the APE of 2.6 KVA. It becomes clear 
that the reactive power compensation using the capacitor bank has reduced the necessary KVA rating of the 
APE. The reduction was approximately 1.8 KVA, from 4.4 KVA to 2.6 KVA, which represents a remarkable 
41% size reduction of the power electronic converter needed for the shunt APE. It is as well important to 
remark that for each particular case, the reduction of the necessary KVA rating in the shunt APE converter 
will depend on the amount of reactive power compensation to be achieved with the capacitor bank during the 
process of transforming the shunt APE compensation scheme into the SHE compensation scheme. The power 
flows during the shunt hybrid filter compensation. Now the capacitor bank is supplying the reactive power 
consumed by the load. As a consequence of a better tracking control over the filtering currents, the losses 
have been significantly reduced. The power losses in the shunt APE were reduced by 62%, i.e. from 184 W 
to 70 W. 


Dsaeca 
Ts= 5e-06s. 


pnvergui 



Eigure 1. Hybrid Shunt Active Power Eilter 
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Figure 2. Non Linear Load 
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Figure 3. Shunt Active Power Filter 
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Under this SHF compensation scheme, the sinusoidal load voltages observed in Figure 3(a) are the 
terminal voltages across the capacitor. Finally, to complete the analysis of results, the currents through the 
capacitor. Where it can be observed that the high frequency ripple current has been drained-out thorough the 
capacitor as it was expected and specified in the control design. The total harmonic distortion of the capacitor 
current (ice) was found to be 3.1%, barely increasing the ram’s current, due to the harmonic ripple current. 


5. CONCLUSION 

An integrated control for a shunt hybrid filter compensation scheme has been proposed. This 
integrated control is based on the design of a LQR-based current controller, for a shunt hybrid filter, which 
has been specifically designed and proposed in this contribution for this SHF scheme of compensation in 
particular. At the end of the process, the structure of the SHF is composed by a shunt Active Power Filter and 
a shunt capacitor bank. Obtained results have shown a remarkable performance on the elimination of 
harmonic distortion (null THD index achieved by the shunt APF) and the compensation of reactive power 
achieved with a conventional shunt capacitor bank is also companied by using the shunt capacitor bank as a 
high-pass ripple filter. Therefore, the harmonic ripple current propagation into the electric system has been 
restricted to a given path... This control action has a very important signification because it provides a sound 
solution to the potential presence of parallel-resonance associated to the shunt active filtering process, as well 
as it avoids the harmonic ripple propagation into the electric network, e.g., to the load voltages at the point of 
common coupling. It has been demonstrated that the proposed control provides a cost-effective solution for 
shunt active power filtering since transforming the shunt APF into a shunt hybrid filter topology reduces the 
overall size and cost of the shunt APF due to the size reduction on the KVA rating needed by the electronic 
power converter of the shunt APF. 
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